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ABSTRACT 



Coronary heart disease (CHD) is a complex and 
heterogeneous cardiovascular disease. There are many 
genome-wide association studies (GWAS) performed 
worldwide to extract the causative genetic factors. 
Moreover, each population may have some excep- 
tional genetic characteristic. Thus, the background 
of our study is from the previous Lithuanian stud- 
ies (the LiVicordia Project), which demonstrated the 
differences of the atherosclerosis process between 
Lithuanian and Swedish male individuals. 

In this study we performed GWAS of 32 families 
of Lithuanian origin in search of significant candidate 
genetic markers [single nucleotide polymorphisms 
(SNPs)] of CHD in this population. After careful 
clinical and biochemical phenotype evaluation, the 
-770K SNPs genotyping (Illumina HumanOmniEx- 
press-12 vl.O array) and familial GWAS analyses 
were performed. 

Twelve SNPs were found to be significantly as- 
sociated with the CHD phenotype (p value <0.0001; 
the power >0.65). The odds ratio (OR) values were 
calculated. Two SNPs (rs 1 7046570 in the RTN4 gene 
and rs 11743737 in the FBXL17 gene) stood out and 
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may prove to be important genetic factors for CHD 
risk. Our results correspond with the findings in other 
studies, and these two SNPs may be the susceptibility 
loci for CHD. 

Keywords: Atherosclerosis; Coronary heart 
disease (CHD); Genome-wide association study(ies) 
(GWAS); Transmission disequilibrium test (TDT). 

INTRODUCTION 



Coronary heart disease (CHD), also called coro- 
nary artery disease, is a complex and heterogeneous 
cardiovascular disease (CVD). It belongs to a group of 
atherosclerotic CVD that is defined as a chronic dis- 
order which develops insidiously throughout life and 
usually progresses to an advanced stage by the time 
symptoms occur [1]. The critical underlying process 
of pathogenesis is atherosclerosis (AS) that, in itself, 
is a multifactorial and peculiar condition. There are 
a number of known controllable and uncontrollable 
factors, one of the last-mentioned is genetic, named 
as strong family history of premature CVD [2]. There 
are many genome-wide association studies (GWAS) 
performed worldwide to determine the main genetic 
factors that could be used for CVD identification and 
creation of useful tests for effective diagnosis, prog- 
nosis and treatment. Regrettably, the genetic factors 
and their importance are not yet sufficiently applied 
in clinical practice [1]. Moreover each population 
may have some exceptional genetic characteristic that 
does not necessarily correspond with results from 
other studies. 
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The background of our study is from the previ- 
ous Linkoping- Vilnius CHD risk assessment study [3], 
which demonstrated the differences of atherosclerotic 
process between Lithuanian and Swedish male individu- 
als. Subsequently, other study aimed to identify potential 
genetic markers associated with AS and CHD in the 
Lithuanian population [4]. The results lacked significant 
values for strong association of single nucleotide poly- 
morphisms (SNPs) and disease. Novel genotyping tech- 
niques and platforms provided an improved opportunity 
for a more precise analysis of whole genome variation 
associated with human complex diseases. Thus, in this 
study we performed the GWAS in 32 families of Lithu- 
anian ethnicity in search of significant genetic mark- 
ers (SNPs) of CHD that may elucidate the underlying 
specificity of AS in this population. 

MATERIALS AND METHODS 



All study protocols were approved by the Vil- 
nius Regional Biomedical Research Ethics Committee 
(No. 158200-ll-255-067LP2;2010-ll-05).Informed 
consent was obtained from all individuals who par- 
ticipated in the study. 

Individuals and Phenotype Definition. Accord- 
ing to the Department of Statistics of Lithuania, in 
2011, the total number of inhabitants was 3,030,200, 
and 17.3% of these were of non Lithuanian ethnic- 
ity. Considering the demographics and the fact that 
trios with CHD patients are generally rare, in total 
we collected 32 relevant families (trios of premature 
CHD patients and their parents) from different regions 
throughout Lithuania, i.e., 96 individuals in total. The 
patient group was represented by 3 1 males and only 
one female. Thus, conclusions from our study are 
suitable to extrapolate only for the male population. 

Patients were clinically examined at the Clinic 
of Cardiology and Angiology, Faculty of Medicine, 
Vilnius University, Vilnius, Lithuania. The patients' 
clinical phenotype was evaluated after assessment 
of anthropometrical measurements, clinical and in- 
strumental examination, and laboratory biochemical 
testing. Information about CHD risk factors, other 
diseases and treatment was obtained during the con- 
ventional anamnesis. 

Patient recruitment criteria was as follows: men 
aged 35-55, women aged 35-65; individuals who were 
experiencing acute coronary syndrome for the first 
time in their lives, who were hospitalized at an inten- 



sive cardiology unit for myocardial infarction (MI) 
with or without Q wave or unstable angina pectoris 
(confirmed by examination of common electrocardio- 
graphic and/or coronographic changes, assessment of 
cardio-specific markers); any previous or current evi- 
dence of significant atherosclerotic CHD (MI, percu- 
taneous coronary angioplasty, coronary artery bypass 
graft or coronary angiography with hemo-dynamically 
significant stenosis). Anatomical vascular changes 
were confirmed by non invasive methods, evaluating 
the presence of the atherosclerotic plaques as well as 
examining the arterial stiffness and endothelial func- 
tion. Patient exclusion criteria was as follows: diabetes 
treated with insulin, kidney function deficiency, III-IV 
functional class of heart deficiency, tumors (except 
skin basalioma), alcoholism and other social factors 
that may influence the study results. 

Metabolic syndrome was diagnosed by determin- 
ing three or more criteria described in the NCEPATP 
///program [5]: waist circumference in men >1 02 
cm, in women >88 cm, triglycerides >1.7 mmol/L, 
high-density lipoprotein cholesterol (HDL-C) <1 .0 in 
men and <1 .2 in women, blood pressure >130/85 mm 
Hg, fasting glucose >5.6 mmol/L or type 2 diabetes. 

Evaluation of biochemical phenotype included 
inflammatory and metabolic markers, participating in 
the patho-genesis of atherosclerosis: C-reactive pro- 
tein (CRP), Hb A |c , lipoprotein [Lp(a)], apolipoprotein 
Al (ApoAl), apolipoprotein B (ApoB), ratio ApoB/ 
Al , lipids, oxidized low-density lipoprotein (oxLDL) 
homocysteine, interleukin-6 (IL-6), fasting glucose in 
plasma, potassium (K), sodium (Na), urea, creatinin. 

Genotyping. The DNA samples were extracted 
from the peripheral venous blood using phenol-chlo- 
roform DNA isolation or TECAN Freedom EVO® 
platform (TECAN Group Ltd., Mannedorf, Switzer- 
land) protocols. All DNA samples were genotyped 
according to the Illumina® provided protocols on 
Illumina HiScanSQ™ genetic analysis platform us- 
ing an Illumina HumanOmni Express-12 vl.O array 
comprised of -770K SNP markers (Illumina Inc., San 
Diego, CA, USA). 

The primary genotyping results were visualized, 
inspected and prepared for further analysis by using 
the GenomeStudio software (Illumina Inc.). Subse- 
quent data quality control, filtering and analysis were 
performed by using the PLINK vl.07 software [6] 
integrated in to the BC|Gene platform (Biocomputing 
Platforms Ltd., Espoo, Finland). 
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Statistical Analysis. The transmission-disequi- 
librium test (TDT) was used to determine the asso- 
ciation between SNPs and CHD phenotype using the 
McNemar test [7]. Single nucleotide polymorphisms 
included in the TDT analysis met the following crite- 
ria: 1) minor allele frequency (MAF) greater than 0.01 
(MAF >0.01); 2) missingness rate smaller than 0.1 
(GENO <0. 1); 3) the p value (O.0005) of the Hardy- 
Weinberg equilibrium test; also if the subjects' missing 
genotype rate was lower than 0.05 (MIND <0.05). For 
the computation of the empirical power, the frequency 
of the informative transmission of disease alleles and 
the number of informative transmissions of marker 
alleles was used. The binomial distribution was used 
for the approximation of the TDT power [8]. The 10 6 
random variables distributed according to the TDT 
statistics under alternative hypothesis were generated. 
Odds ratio (OR) with 95% confidence interval (CI 
95%) was calculated for selected SNPs. The level of 
significance was set at a = 1 0 -4 . For multiple testing, 
the correction adaptive permutation procedure was 
used. Maximum number of permutations was 1 0 6 . The 
statistical analysis was performed with the program 
R (v2. 15.3) (R Foundation for Statistical Computing, 
Vienna, Austria). 



RESULTS 



According to the Illumina Inc. protocol guide- 
lines, all of the samples except one, were of good 
quality and had been properly processed (call rate 
>98; LogRDev <0.3; coincidental sex list file created). 
At the beginning of the analysis there were 731,412 
SNPs genotyped in the group of 96 individuals. Af- 
ter the data filtering procedure, two individuals were 
removed from further analysis for low genotyping 
(MIND >0.05); 25,293 heterozygous genotypes were 
excluded from analysis because the second allele of the 
genotype was missing; 298 SNPs were excluded based 
on the Hardy-Weinberg equilibrium test (p >0.0005); 
2528 SNPs failed missingness test (GENO >0.1); 
82,552 SNPs failed frequency test (MAF <0.01); 591 
SNP were not used because of homogeneity over all 
individuals. After the final frequency and genotyping 
pruning, 646,445 SNPs in 3 1 patient and 63 parents 
were included for further association analysis. 

Twelve SNPs were found to be significantly as- 
sociated with CHD phenotype with p values smaller 
than 0.0001 . The SNPs annotation (transmitted allele, 
chromosomal position, gene, gene function) along 
with the x 2 , P value, OR and empirical power based 
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Figure 1. Manhattan plot. Single nucleotide polymorphism distribution with logarithmic transformation of the 
p value. The horizontal line depicts the selected significance level (a = 0.0001). The x axis represents the SNPs' posi- 
tions according to each chromosome, the y axis shows -loglO(P) values of the SNPs. 
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Table 1. Statistically significant single nucleotide polymorphisms associated with coronary hearth disease. 



Chromosome 


SNP 


Transmitted 
Allele 


Position 


Gene 


Function 


X 2 


p Value 


OR (95% CI) 


Power 


2pl6.1 
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1 Q 1A 
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U.OZ 


jqjl.l 
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1 j4j / jyyl 


LOC340073; 
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U. Ij (U.U4-U.4Z) 


A £Q 
U.00 


f «") 7 3 

3qzl.3 


mi* 7 7 7,7 2 72 7 

rsii 743 73 / 


71 


IV /ZojjoV 


I7DVT 7 7 

rlSALl 7 


SCF complex; protein 
ubiquitination 


If ID 

13.39 


0.77x10 


T/7 /I 2/1 f?i 

/.O/ fz.3U-z3.33) 


fi AZ 
U.OJ 


9q22.23 


rsl 08 19695 


G 


98762901 






15.70 


7.439xl0- 5 


0.10(0.02-0.41 


0.69 


12pll.23 


rsll048567 


A 


26689621 


ITPR2 


intracellular Ca 2+ 
relaxation 


17.19 


3.38xl0 5 


0.05 (0.01-0.37 


0.81 


20qll.21 


rs6141273 


A 


29904377 






16.20 


5.70xl0 5 


0.05 (0.01-0.39) 


0.77 


20pl2 


rsl321936 


G 


12826999 






16.03 


6.23x10 s 


0.18(0.07-0.46 


0.64 




Possible Artefacts 


lq32.1 


rsl2734338 


G 


202469723 


PPP1R12B 


regulation myosin 
phos-phatase activity; 
aug-ments Ca 2+ 
sensitivity of the 
contractile apparatus 


26.13 


3.19x10-' 


29.0 (3.95-212.9) 


0.997 


4q25 


rs3853444 


G 


111734136 






20.83 


5.023x10-" 


0.13 (0.05-0.37) 


0.880 


15q25.2 


rs3883013 


G 


85088657 


UBE2Q2P1 


psuedogene 


26.13 


3.19x10-' 


29.0 (3.95-212.9) 


0.997 



Ca 2+ : calcium ion. Out of 12 significant SNPs, at least two (italicized and bold rows) had promising OR values in addition to power and p values. 



on the sample size calculations, are presented in Table 
1 . The SNPs are annotated according to the National 
Center for Biotechnology Information (NCBI) dbSNP 
and Gene databases [9]. 

The acceptable power values were greater than 
or equal to 0.65, and thus fell partly into the desired 
range between 0.8 and 0.95 [10]. Only the power 
value of the significant SNP rsl 32 1936 diverged and 
was excluded from further evaluation. 

The OR values in Table 1 show the size of the 
effect. The greater the deviation of OR is from the 
value of 1, the more significant the test is. 



As can be seen from the Manhattan plot (Figure 
1), there are three significant markers (rs 12734338, 
rs3 8830 1 3, rs3 853444) that do not have the significant 
adjacent SNPs (according to the nucleotide's posi- 
tion), i.e., the correlation is absent. Thus, these mark- 
ers could be artefacts. We could also suspect that not 
all of the adjacent SNPs were genotyped. This is more 
likely to happen with rs3853444, as there were two 
adjacent SNPs that were excluded from the analysis. 

It was previously mentioned that the study group 
included only male patients and their parents. A male 
could possess the SNP allele on either his X or Y 
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chromosome and this affects the analysis algorithm. 
Each transmission from a heterozygous mother to a 
male offspring should be given twice the weight of 
a transmission to a female offspring [11]. Thus, the 
standard TDT appears to be unsuitable for the analy- 
sis of SNPs in sex chromosomes and eventually sex 
chromosomes were excluded from the analysis. 

DISCUSSION 



Our aim was to identify the particular genetic 
factors for Lihuanian CHD patients. In order to avoid 
population stratification we planned the familial 
GWAS and subsequently performed TDT analysis. 
Out of 12 significant SNPs, at least two (italicized 
and bold rows in Table 1) had promising OR values 
in addition to power and p values. Despite relatively 
wide OR intervals, caused by a modest sample size, 
these SNPs may indicate the potential genes that could 
be involved in CHD pathogenesis. 

The SNP rs 17046570 is located in an intron of 
the reticulon 4 coding gene (RTN4) on chromosome 
2. Retic-ulons are associated with endoplasmic re- 
ticulum, and are involved in neuroendocrine secretion 
or membrane trafficking and apoptotic processes. In 
particular, reticulon 4 has been identified as a potential 
inhibitor of central nervous system regeneration by 
means of the inhibition of neuron outgrowth [11,12]. 
Common RTN4 variants that are associated with 
schizophrenia in the Japanese population [13] and 
also blood lipid phenotypes [14], are cited. RTN4 is a 
candidate gene associated with vascular cell apoptosis 
and AS modulation [12]. It is thought to participate in 
vascular remodeling and is a considerable new factor 
for atherogenesis process [15-17]. It was also stated 
that reticulons may be factors that mediate between 
the apoptosis and AS processes [12]. Thus, our results 
are consistent with these findings. 

Another SNP on chromosome 5, rsl 1743737, 
is located in an intron of the F-box and leucine-rich 
repeat protein 1 7 coding gene FBXL1 7. The FBXL 1 7 
protein has an F-box that is a 40 amino acid motif 
typical for F-box containing proteins. F-box contain- 
ing proteins together with culin and SKP1 (S-phase 
kinase-associated protein 1) make up the SCF com- 
plex (SKP1, cullin, F-box containing complex) that 
is a protein ubiquitine ligase [18]. The SCF is a key 
complex in the ubiquitine-proteosome system (UPS) 
that is involved in 70.0-90.0% of protein degradation 
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processes including the degradation of a number of 
proteins important for the cardiovascular system. The 
UPS is also important in the regulation of endothelial 
cell cycle. The effect of oxidative stress on the SCF 
complex may disrupt the function of UPS and in turn 
the function of the endothelium that is regulated the 
by UPS [19]. According to the NCBI Gene database 
review of association results [from National Human 
Genome Research Institute (NHGRI) Catalogue and 
association results submitted to the database of Geno- 
types and Phenotypes (dbGaP)] there are many SNPs 
in the FBXL 17 gene region associated with the vari- 
ous phenotypes including cholesterol, high-density 
lipoproteins, body mass index. These findings do not 
compromise our findings either. Moreover, these sum- 
marized results might show us the complexity and 
universality of the FBXL 17 protein function in the 
pathogenesis of different diseases. It is possible that 
other SNPs that are in linkage disequilibrium with the 
identified CHD associated SNPs were not identified 
during this TDT analysis but may also be involved in 
the development of the disease. 

Conclusions. Our results suggest that the RTN4 
and FBXL1 7 genes may be the susceptibility loci for 
the CHD in the Lithuanian male population. In ad- 
dition, the genotypes of the significantly associated 
SNPs with OR >7 may prove to be informative and 
specific for the genetic risk of CHD evaluation in 
the Lithuanian population and could be taken under 
consideration in further hypothesis validation. 
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